\j 


t  ■n 


I 

C 


SEL-63-011 


Some  Backscatter,  Sweep-Frequency  Observations 

402  964 

by 

S.  N.  Denno*  and  0.  G.  Villard,  Jr. 


c 

i- 

c 

c 


February  1963 

' 


Technical  Report  No.  72 


Prepared  under 

Office  of  Naval  Research  Contract 

Nonr-225(64),  NR  088-003,  and 

Advanced  Research  Projects  Agency  ARPA  Order  196-62 


r  *  -  T  1  A 

nr-? 


*  Supported  by  National  Research  Council  Fellowship 
National  Academy  of  Sciences 


t » ..  ■ 

it.'  '  y 


j  j  ,1 


n 

m 

■  i : 


l  -v;:  a 


RRDIOSIIEnCE  LHB0RRT0RV 

SinniORD  ELGCTROniCS  LRB0RRT0RIES  5 

STRnFORD  URIVERSITV  •  STflnFORD,  criiforrir  > 


ASTIA  AVAILABILITY  NOTICK 

Qualified  requesters  may  obtain  copies  of  this 
report  from  ASTIA.  Foreign  announcement  and 
dissemination  of  this  report  by  ASTIA  is  limited. 


S.  N.  Deruno  and  0.  G.  Villard,  Jr 


7? 

TR  No. 

SBL-63-OU 

February  1963 


ERRATA 


Page 

vi 

Location 

Should  Read 

Line  3  from 
bottom 

at  the  assumed  position  of  the  irregulatity 
...53 

6 

Line  3  from 
bottom 

|l/(dP'/d0o)| .  It  is  obvious  that... 

11 

Line  5  below 

Eq.  (l6) 

path.  This  function  provides  a  rough  approxl 
mation. . . 

Line  2  below 

Eq.  (17) 

with  this  angle  and  occurs  at  an  additional 
normalized  height  of 

lit 

Line  6  below 

Eq.  (19) 

transmitter,  viz.,  0°,  +30°,  +60°,  and  -6o° 
for  b  equal  to. . . 

2  It 

Line  8  from 
bottom 

...of  distance  traveled  by  each  ray  on  echo 
intensity,  the  num- 

30 

Line  12  from 
bottom 

... in  the  region  of  maximum  frequency,  no 
longer  bears . . . 

33 

Fig.  21  A-l 
above  top 
strip 

22  June  1962 

1*5 

Line  4 

frequency  range  of  15  -  24  Mc/s,  show,  on 
the ... 

53 

Last  line  of 
legend 

ASSUMED  POSITION  OF  THE  IRREGULARITY. 

63 

Eq.  (D.3) 

Eq.  (D.6) 

BC  =  PC  cos  0. 

...sin2  0  =  -2  FA  ’  PC  cos  a... 

SEL- 63-01 1 


SOME  BACKSCATTER, 


SWEEP- FREQUENCY  OBSERVATIONS 


by 


S.  N.  Denno*  and  0.  G.  Villard,  Jr. 


February  1963 


Reproduction  in  whole  or  in  part¬ 
is  permitted  for  and  purpose  of 
the  United  States  Government. 


Technical  Report  No.  72 

Prepared  under 

Office  of  Naval  Research  Contract 
Non  r- 225(64),  NR  088-003,  and 
Advanced  Research  Projects  Agency  ARPA  Order  196-62 

Supported  by 

National  Research  Council  Fellowship 
National  Academy  of  Sciences 


Radioscience  Laboratory 
Stanford  Electronics  Laboratories 
Stanford  University  Stanford,  California 


ABSTRACT 


Analysis  of  the  sweep-frequency,  ground -hackscatter  echo  is  attempted 
by  the  application  of  ray  theory.  The  variation  of  the  time  delay  of  the 
various  portions  of  the  radar  echo  is  studied  through  two  different 
methods:  l)  calculation  of  the  equivalent  group  path,  and  2)  tracing 
the  actual  group  path  by  means  of  an  analog  computer.  Special  attention 
is  given  the  two  limiting  rays --the  tangent  ray  and  the  minimum-time-delay 
ray — in  the  study  of  their  behavior  with  frequency,  since  they  determine 
the  shape  of  the  sweep-frequency  pattern.  Both  approaches  agree  well  with 
experimental  observations.  The  effects' of  the  layer  parameters  on  the 
shape  are  briefly  investigated;  these  include  the  layer  model,  its  thick- 
ness,  and  the  presence  of  horizontal  gradients. 

A  study  of  the  striations  that  appear  within  the  Sweep-Frequency 
Echo  (SFE)  is  also  reported.  A  description  is  given  of  continuous  sound¬ 
ings  made  over  a  limited  period  in  the  summer  of  1962  in  the  eastward 
and  westward  directions  from  Stanford,  California,  for  the  purpose  of 
determining  their  time  and  frequency  characteristics.  Two  distinct  types 
are  separated  according  to  the  variation  of  their  time  delay  with  fre¬ 
quency.  With  the  aid  of  the  above  analysis,  it  is  shown  that  the 
frequency-dependent  striations  may  be  caused  by  layer  stratification  or 
by  the  presence  of  not  too  intense  electron-density  "blobs"  embedded  in 
the  F -layer,  while  those  that  are  frequency-independent  are  the  result 
of  highly  reflective  irregularities  in  the  propagation  path.  The  loca¬ 
tion  of  these  irregularities  is,  so  far,  a  matter  of  speculation,  but 
it  is  thought  that  the  E-layer  is  a  quite  possible  one,  and  that  some  of 
these  irregularities  may  be  identified  with  ordinary  sporadic  E  and  also 
with  the  middle-latitude  variety  of  field-aligned  irregularities.  The 
geometrical  conditions  are  explored  under  which  the  latter  kind  may  be 
illuminated  by  an  F-layer-propagated  ray  in  the  downward  part  of  its 
path  while  satisfying  normal  incidence  with  their  long  axis.  Further 
work  is  reconsnended  to  pin  down  the  location  of  these  irregularities  and 
to  extend  knowledge  of  their  characteristics. 
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ground  distance  between  transmitter  and  subsolar  point 
velocity  of  light  in  free  space 

maximum  frequency  at  which  ray  will  be  reflected  at  peak 
electron  density  in  layer  (critical  frequency) 

wave  frequency  at  oblique  propagation 

wave  frequency  at  vertical  propagation 

true  height  of  ray  reflection  in  the  F  layer  (in  ray  tracing 
h  and  x  were  used  as  the  rectangular  coordinates  of  the  ray 
path) 

height  of  bottom  of  ionosphere  (F  layer  of  parabolic  distribu¬ 
tion) 

virtual  height  of  ray  reflection 

correction  factor  for  curvature  of  ionosphere 

group  velocity  of  propagated  wave 

ground  distance  of  ray  from  subsolar  point 

half  thickness  of  ionosphere  (F  layer  of  parabolic  distribu¬ 
tion) 

normalized  height  used  in  Chapman  formula  for  ion-density 
distribution,  equal  to  (h  -  h^nv)/H 

angle  between  ray-path  plane  and  magnetic -north  direction  at 
a  point  in  the  path  at  E-layer  height 

geomagnetic-declination  angle  at  E-layer  height.  Also  ground 
distance  covered  by  ray  between  transmitter  and  scattering 
point  on  ground 

vertical-distribution  function  of  electron  density  in  F  layer 
geomagnetic-dip  angle  at  E-layer  height 
electron-number  density  in  ionosphere  (F  layer) 
peak  electron-number  density  in  F  layer 

total  length  of  equivalent  group  path  of  ray  between  transmitter 
and  scattering  point  on  ground,  i.e.,  oblique  distance  of 
scattering  point 


SYMBOLS  (Cont'd) 


P'  length  of  oblique  distance  corresponding  to  minimum  delay  time 

mllx  (close  to  skip  distance) 

R  earth's  radius  (6367  km) 

S  total  length  of  true  group  path  of  ray  between  transmitter  and 

scattering  point  on  ground 

S(x)  horizontal-distribution  function  of  electron  density  in  F 
layer  along  ray  path 

S(y)  horizontal-distribution  function  of  electron  density  in  F 
layer  in  a  meridian  plane  containing  subsolar  point  and  a 
point  in  ray  path 

T  total  one-way  group  delay  of  ray  between  transmitter  and  scatter¬ 

ing  point  on  ground 

a  angle  between  the  two  meridian  planes  containing  propagated 

ray  and  a  geomagnetic  field  line  at  E-layer  height,  respectively 

initial  vertical  angle  of  propagated  ray 

p  vertical  angle  of  downcoming  ray  at  E-layer  height  after 

p  reflection  by  F  layer 

p.  phase-refractive  index  at  any  point  in  ionospheric  layer 

supporting  propagation 

p  normalized  frequency,  equal  to  f/f  ,  where  f  is  frequency  of 

propagated  wave.  For  ground  backs cattering  f  *  f  . 

ob 

0  angle  of  incidence  of  ray  at  the  point  of  virtual  reflection  in 

0  the  ionosphere 

X  zenith  angle  of  sun's  rays  incident  on  ionosphere 

0  angle  at  earth's  center  subtended  by  arc  distance  D 
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I.  INTRODUCTION 


The  work  reported  here  is  an  analysis  of  the  ground -backs cat ter  echo, 
with  particular  reference  to  the  effect  of  continuous  variation  of  fre¬ 
quency  and  the  appearance  of  certain  echoes  that  seem  to  be  unaffected  by 
this  variation,  thus  appearing  as  constant-range  echoes. 

It  is  now  well  established  that,  when  a  pulse  of  radio  frequency  is 
propagated  obliquely  via  the  ionosphere,  an  echo  is  received  after  a 
certain  time  delay  as  a  result  of  scattering  from  irregularities  on  the 
distant  earth's  surface,  be  it  land  or  sea,  providing  that  one  or  more 
of  the  ionospheric  layers  can  support  such  propagation.  If  the  pulse 
length  is  of  the  order  of  1  msec,  the  nature  of  the  scattering  irregu¬ 
larities  seems  to  have  a  secondary  effect  on  the  strength  of  the  pattern 
formed  by  the  received  echoes  [Ref.  1],  Thus,  whether  the  beam  is  directed 
towards  water  or  dry  land,  the  pattern  is  governed  primarily  by  the  pre¬ 
vailing  structure  of  the  ionosphere.  Furthermore,  because  of  the  random 
distribution  of  the  scattering  irregularities,  the  received  echo  is 
phase  incoherent j  but  nevertheless,  it  is  possible  to  employ  the  princi¬ 
ples  of  simple  ray  theory  in  the  analysis  of  the  transmission  mechanism. 
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II.  CHARACTERISTICS  OF  THE  BACKSCATTER  ECHO 


A.  THE  REFLECTED  ECHO  AND  ITS  PRESENTATION 

When  a  radar  pulse  is  transmitted  from  an  isotropic  radiator,  the 
shape  of  the  returning  echo,  after  being  reflected  twice  by  the  iono¬ 
sphere  and  scattered  once  by  the  ground,  is  governed  by  the  time  delay 
of  each  ray  component  and  the  density  of  these  components  at  any  particu¬ 
lar  moment.  The  original  shape  of  the  pulse  is  thus  completely  lost 
[Ref.  2].  The  analysis  follows  each  ray  separately  along  its  group  path, 
treating  it  as  an  individual  pulse  of  infinitesimal  duration  with  appli¬ 
cation  of  the  laws  of  oblique  transmission.  Hence  the  characteristics 
of  the  ionosphere  present  along  the  ray  path  have  a  primary  effect  on 
this  shape. 

The  backscatter  echo  may  be  presented  in  an  A-scope  display  relating 
received  field  strength  and  time  delays,  or  by  using  the  electron-beam 
modulation  of  a  cathode-ray  tube  as  a  third  dimension.  The  above  rela¬ 
tionship  can  be  extended  to  show  the  variation  of  the  echo  as  a  function 
of  time  or  frequency. 

B.  GROUP  PATH  AND  TIME  DELAY 

The  aforementioned  display  gives  the  delay  of  the  echo,  2T,  which  is 
the  time  taken  for  the  pulse  to  make  the  round  trip  before  it  is  received 
as  an  echo.  It  is  assumed  here  that  the  actual  group  path  S  followed 
by  the  ray  in  the  outgoing  direction  is  the  same  as  that  followed  in 
the  return  direction;  hence,  for  simplicity,  we  shall  consider  in  our 
analysis  one-way  conditions  only.  Since  the  ionosphere  is  geographically 
non-homogeneous,  the  phase-refractive  index  p  and  therefore  the  velocity 
of  propagation  will  be  functions  of  space.  For  a  group  velocity  v 

6 

the  time  delay  is  given  by 


and,  when  the  effects  of  the  earth's  magnetic  field  and  collisional 
damping  are  neglected,  the  time  delay  can  be  expressed  in  terms  of  the 
phase  refractive  index  as 
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where  e  is  the  velocity  of  electromagnetic  waves  in  vacuum. 

The  concept  of  the  equivalent  path  is  very  useful  in  sounding  work 
and,  in  this  case,  the  equivalent-group  path  P1  simply  replaces  the 
integral  term  in  Eq.  (2).  Thus, 


Figure  1  shows  the  geometry  of  the  equivalent  path  for  curved  earth 
condition.  B  is  a  point  of  assumed  specular  reflection  at  an  equiva¬ 
lent  height  h'  in  the  ionosphere.  Under  flat-earth  conditions,  the 
equivalence  theorem  relates  P'  and  h'  by  the  equation 


h* 


-  f  CO,  0o, 


(■*) 


FIG.  I.  GEOMETRY  OF  EQUIVALENT  RAY  PATH. 


ionosphere  and  is  also  the  half-vertex  angle  of  the  equivalent  triangular 
path. 

The  curvature  of  the  earth  and  ionosphere  requires  modification  of 
the  plane  geometry  as  well  as  the  so-called  secant  law  on  which  the 
equivalence  theorem  is  based.  Smith  [Ref.  3]  has  shown  that,  in  the  case 
of  the  ionosphere  following  the  curvature  of  the  earth's  surface  (radius 
R)  with  a  lower  boundary  at  a  constant  height  hQ  above  it,  the  true 
relationship  between  the  oblique-incidence  frequency  f  ^  and  the 
equivalent  vertical-incidence  frequency  f  for  which  h'  is  the 
vertical  height  is  given  by 

fv  =  fob  COS  0o  [L  -  (RhI  hQ))tan2  *o]  2  (5a> 

=  f ob  (cOS  V  k0o>  (5b) 

where  h  is  the  real  height  at  which  the  actual  ray  trajectory  reaches 
at  its  apex. 

Note  that  the  angle  0q  remains  to  be  half  the  vertex  angle  of  the 
triangle  ABC,  but  in  this  case  it  is  smaller  than  the  angle  of  incidence 
at  the  bottom  of  the  ionosphere. 

It  will  also  be  noticed  that  the  effect  of  the  factor  k^Q  in  Eq. 
(5b)  is  to  lower  the  equivalent  vertical-incidence  frequency  or  what 
amounts  to  the  lowering  of  the  virtual  height  of  vertical  incidence, 
which  can  readily  be  seen  from  the  curved  geometry. 

It  is  sometimes  more  convenient  to  express  the  above  relationships 
in  terms  of  the  take-off  angle  for  the  curved-earth  condition,  as 

shown  in  Appendix  A,  the  equivalent  path  is  given  by 

P'  =  2{[(R  +  h')2  -  R2  cos2  Pq]5  -  R  sin  3  ).  (6) 

For  a  layer  of  parabolic  electron  distribution  of  half  thickness  y  , 

m 

the  value  of  h'  can  be  derived  from  Eq.  (A.  12)  of  Appendix  A,  which 
is  in  a  form  amenable  to  simple  solution  by  a  digital  computer. 
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C.  THE  EQUIVALENT  PATH  APPROACH 


The  study  of  the  characteristics  of  the  ground  backs catter  echo  Is 
first  made  by  utilizing  the  concept  of  the  equivalent-group  path,  for 
which  it  is  necessary  to  assume  a  certain  ion-distribution  model  of  the 
reflecting  layer. 

A  useful  graphical  method  for  this  study  is  that  developed  by  Peter¬ 
son  [Ref.  4],  which  shows,  in  combined  rectangular-polar  coordinate 
systems,  the  relationship  between  take-off  angle  equivalent -group 
distance  or  oblique  path  P1,  virtual  height  of  reflection  h',  and 

horizontal  ground  distance  to  the  point  of  ground  reflection  D.  The 

* 

curve  relating  these  quantities,  which  may  be  called  the  reflectrix, 
is  drawn  parametric  with  the  operating  frequency  f  ^  which  is  usually 
normalized  with  respect  to  the  vertical-incidence  critical  frequency  of 
the  ionospheric  layer  f  ^ .  This  treatment,  however,  was  confined  to  the 
case  of  flat  earth  and  flat  ionosphere,  and  its  accuracy  is  therefore 
limited  to  short  ranges.  For  longer  ranges,  the  curved  geometry  shown 
in  Fig.  1  should  be  taken  into  consideration,  and  Peterson's  construction 
can  still  be  applied  with  the  same  coordinate  system,  for  a  given  iono¬ 
spheric  electron  distribution.  The  shape  of  the  reflectrix,  i.e.,  the 
locus  of  the  extremity  of  the  line,  (in  this  case  curved)  representing 
the  length  of  the  equivalent  path,  however,  will  be  quite  different. 

Alternatively,  the  construction  can  be  made  more  realistic  by  using 
a  circular  coordinate  system  as  shown  in  Fig.  2.  The  ground  range  D 
corresponding  to  any  point  on  the  reflectrix  will  be  measured  by  the 
radial  projection  along  a  circle  of  radius  2R,  where  R  is  the  earth's 
radius,  thus  keeping  the  angular  relationship  of  the  transmission  geom¬ 
etry  the  same,  with  a  linear  representation  of  distances  h'  and  P'. 
Good  approximation  may  be  obtained  by  simply  drawing  the  plane-earth 
data  on  the  circular  coordinates. 

Since  the  length  of  the  equivalent  path  P’  is  proportional  to 

the  time  delay  T  (Eg,  1)  the  construction  in  Fig.  2  facilitates  the 
* 

The  term,  originally  devised  by  Le jay  and  Lepechinsky  [Ref.  5],  applied 
to  the  locus  of  the  virtual  point  of  the  first  reflection  in  the  iono¬ 
sphere.  As  it  will  be  used  here,  the  reflectrix  involves  double  the 
distances  and  the  same  angles. 
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FIG.  2.  GRAPHICAL  CONSTRUCTION  SHOWING  THE  VARIATION 
OF  THE  EQUIVALENT  GROUP  PATH  P'  OF  A  RAY  WITH  THE 
ELEVATION  ANGLE  AT  DIFFERENT  NORMALIZED  FRE¬ 

QUENCIES  p. 


study  of  the  shape  of  the  backscatter  echo  by  following  the  variation 
of  P1  along  any  of  the  reflectrices  shown.  For  any  normalized  fre¬ 
quency  P,  the  width  of  the  echo,  i.e.,  its  duration,  will  fall  between 
points  such  as  N  and  M.  The  former,  at  the  tip  of  the  reflectrix, 
is  the  point  of  minimum  time  delay,  and  the  latter  corresponds  to  the 
maximum  length  of  path  at  zero-angle  propagation.  Furthermore,  the  rate 
of  change  of  the  length  of  P1  with  respect  to  the  angle  (3q  gives  a 
measure  of  the  amplitude  of  the  echo,  so  that,  at  any  value  P1  or  its 
equivalent  delay  T,  the  echo  intensity  depends  to  a  first  order  on 
l/(dP'/d£o).  It  is  obvious  that  the  point  of  minimum  time  delay  produces 
tne  highest  amplitude  in  the  echo,  resulting  from  what  is  known  as 
minimum-time -delay  focussing. 
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It  will  also  be  noticed  that,  unlike  the  flat  earth  case,  the 
intensity  of  the  echo  corresponding  to  zero  take-off  angle  has  a  finite 
value,  as  can  be  seen  from  the  value  of  dP'/d£0  (derived  in  Appendix  A), 
which  reaches  a  limiting  value  of  (-2R)  because  of  the  low  penetration 
at  such  low  angles  (Eq.  A. 15).  This  will  later  be  referred  to  as  the 
limit,  or  tangent-ray  echo. 

Using  this  construction,  Eq.  (6)  is  plotted  for  a  typical  parabolic 
layer  of  maximum  ionization  height  of  300  km  and  semi-thickness  of  100 
km.  The  result  is  shown  in  Fig.  3-  No  correction  was  made  in  this 
case  for  the  curvature  of  the  ionosphere  at  the  point  of  assumed  re¬ 
flection.  Thus  the  factor  was,  for  simplicity,  made  equal  to 

unity. 

In  the  type  of  display  used  in  the  present  experiments --namely, 
range  vs  frequency,  with  intensity  modulation  representing  the  strength 
of  the  received  echo — the  pattern  produced  (called  the  Sweep  Frequency 
Echo,  SEE)  for  a  given  frequency  range  is  bounded  by  two  extreme  echoes 
identifying  the  two  rays  referred  to  earlier- -namely,  the  minimum-time- 
delay  and  the  tangent  ray.  It  is  therefore  useful  to  single  out  these 
rays  and  show  their  variation  with  frequency,  which  is  presented  in  the 
plot  of  Fig.  4,  taken  from  the  construction  of  Fig.  3  after  suitable 
conversion  of  the  equivalent -group  paths  into  one-way  time  delays. 

D.  THE  RAY-TRACING  APPROACH 

1.  General 

A  more  realistic,  though  involved,  approach  to  the  determination 
of  time  delay  and  the  prediction  of  the  shape  of  the  received  echo  is  by 
means  of  ray- tracing  techniques.  The  actual  path  of  the  electromagnetic- 
wave  energy,  which  is  the  true  group  path,  is  followed  along  the  non- 
uniform  electron  concentration  of  the  ionosphere  from  the  point  of 
departure  at  the  transmitter.  Thus  it  is  possible  for  a  study  to  be 
made  when  horizontal  gradients  or,  in  fact,  any  form  of  electron-density 
distribution  exist. 
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USING  EQS.  (7)  AND  (12)  OF  APPENDIX 


FIG.  4.  SHAPE  OF  SWEEP- FREQUENCY ,  GROUND- BACKSCATTERED  ECHO, 
BOUNDED  BY  ITS  TWO  EXTREME  TIME  DELAYS  AS  FUNCTIONS  OF  THE 
NORMALIZED  FREQUENCY.  Parabolic  layer  h0  =  200  km,  ym  =  100  km. 


2.  Ray  Tracings  with  a  ft -Chapman  Layer 

An  analog  computer  was  programmed  to  solve  the  following  equation 
from  ray  theory  [Ref.  6 ],  based  on  Snell's  law  and  Fermat’s  principle  of 
stationarity,  for  curved  earth  and  curved  isotropic  ionosphere  in  a  two- 
dimensional  coordinate  system 


where:  h  =  real  height  above  the  earth's  surface  at  any  arc  distance 
x  from  the  transmitter 

(i  =  phase -refractive  index. 
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In  the  programming  of  the  analog  computer,  two  approximations 
were  made  to  simplify  the  above  equation  for  more  efficient  use  of  the 
computer. 

1.  Since  h  «  R, 

(i  +|)  -l  +  {p,  (8) 

and,  in  the  second  term  of  Eq.  (j), 

(R  +  h)  »  R.  (9) 

2.  Because  of  the  small  variation  of  the  refractive  index  in  the 
horizontal  direction,  the  last  term  of  Eq.  (7)  is  neglected, 
leaving  the  terms  for  the  vertical  gradient  and  the  earth's 
curvature.  However,  the  effect  of  the  horizontal  gradient  of  the 
electron  density  is  not  completely  dropped,  since  it  appears  in 
the  following  evaluation  of  the  refractive  index  in  terms  of 
frequency, 

S-i — "-fe-sUr .  do) 

1.24  x  10  f\ 

OD 

where  N(h,  x)  is  the  electron-number  density  per  cubic  centimeter 
along  any  particular  ray  path. 

If  Nq  is  the  maximum  electron  density  present  in  the  layer,  the 
electron  density  at  any  other  point  (h,  x)  can  be  expressed  in  terms 
of  a  vertical  and  horizontal  components,  F(h)  and  S(x)  respectively, 
of  the  spatial  distribution  function.  Thus, 


Hence, 


N(h,  x)  =  N  •  F(h)  •  S(x)  . 


(11) 


H2  =  1 - ^-3 - 5-  F(h)  S(x)  .  (12) 

1.24  x  10  x  f  . 

ob 


The  concentration  N  will  correspond  to  a  critical  frequency 

0  4  ? 

fQ  Mc/s  at  vertical  incidence  such  that  Nq  =  1.24  x  10  f  ,  giving, 
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(13) 


for  the  refractive  index, 

»i2  -  1  -  F(b)  '  S(x) 

or 

.  i .  (n) 

P 

The  functions  F(h)  and  S(x)  were  produced  in  the  analog 
computer  by  means  of  function  generators  of  the  diode  type.  For  the 
vertical  direction,  a  £ -Chapman  distribution  was  assumed  for  the  F 
layer  [Ref.  7]  with  a  height  of  maximum  ionization  of  300  km  and  a  scale  - 
height  of  100  km.  Thus, 

F(h)  =  exp  [1  -  z  -  exp(-z)],  (15) 

where 

z  -  (h  -  300)/l00.  (l6) 

The  horizontal  gradient  S(x)  of  the  spatial-distribution  function 
is  taken  as  unity  for  the  sample  ray-path  tracings  shown  in  Fig.  5>  A 
cosine  function  was  also  introduced  to  simulate  the  effect  of  the  sun's 
zenith  angle  x  over  'the  various  portions  of  the  ionosphere  along  the 
path.  This  function  provide  a  rough  approximation  to  what,  even  theo¬ 
retically,  is  considered  to  happen  since,  although  the  peak  electron 
density  of  the  P-Chapman  layer  is  proportional  to  the  cosine  of  the 
zenith  angle,  i . e . , 

"o  -  IM.0  •  C0"X'  <17> 

the  height  of  this  peak  does  not  remain  constant  but  rather  increases 
with  this  angle  and  occurs  at 

z  =  In  secx,  -  (17a) 

At 

which  is  ignored  in  the  present  computation. 
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Hence,  at  any  point  along  the  ray  path  whose  projection  on  the 
ground  is  distant  y  km  from  the  subsolar  point,  the  value  of  the 
electron  density  gradient  formed  is 

S(y)  =  cos  (y/R).  (l8) 


In  Appendix  B  the  trigonometry  of  this  general  case  is  given. 
It  is  also  shown  that,  for  the  case  in  which  the  transmitter,  the  sub¬ 
solar  point,  and  the  ray  path  fall  in  the  same  meridian  plane, 


S(x)  =  cos  (y/R)  =  cos  (g  -  | )  ,  (19) 


where  b  is  the  ground  distance  between  the  transmitter  and  the  subsolar 
point.  Note  that  the  frequency  represented  by  p  remains  normalized 
with  respect  to  the  critical  frequency  of  the  maximum  ion  density, 
which  in  this  case  will  be  over  the  subsolar  point.  Ray  tracings  for 
this  case  were  made  with  the  sun  at  different  zenith  angles  from  the 
transmitter,  viz.,  0°,  +30°,  +60°  for  b  equal  to  0,  3333,  6667,  and 
-6667  km,  respectively,  and  are  shown  in  Pigs.  6  through  9.  These 
would  correspond  to  conditions  during  different  hours  of  the  day, 
depending  upon  latitude  and  direction  of  propagation. 

The  ray  tracings  obtained  are  for  different  normalized  frequen¬ 
cies  p  and  initial  angles  of  propagation  given  by 
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b.  From  0°  through  Critical  at  26  ;  P  ~  2.0. 
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c.  From  0°  through  Critical  at  17”;  p  -  2.5. 


d.  From  0°  through  Critical  at  10°;  p  =  3.0. 

FIG.  6.  TRACINGS  OF  RAYS  OF  VARIOUS  TAKE-OFF  ANGLES 
IN  A  /3-CHAPMAN  LAYER.  S(x)  =  cos  [(x-b)/R];  b  =  0  km. 
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.  From  0°  through  Critical  at  51  \  p  -  1.0. 
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b.  From  0°  through  Critical'at  37  \  P  ~  1.5. 
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b.  Between  0  and  21  (Critical);  p 
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c.  Between  0°  and  14°  (Critical);  p  -  1.75. 


FIG.  9.  TRACINGS  OF  RAYS  OF  VARIOUS  TAKE-OFF  ANGLES 
IN  A  /3-CHAPMAN  LAYER  AT  1-deg  INTERVALS.  S(x)  = 
cos  [ ( x—  b ) /Rl ;  b  =  -6667  km. 
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3*  Derivation  of  Time  Delay  and  Its  Relation  to  the  Backscatter  Echo 


The  determination  of  the  group  time  delay  is  obtained  from  the 
computer  loop  by  tapping  off  the  variables  (dh/dx)  and  p.  and  sub¬ 
programming  Eq.  (2l)  (see  Appendix  C) 

h=0  L  /  dh 

,.i  f  L1  *  (g 

C  J  P 

x="> 


> 


1/2 


A  ,, 


/oi  \ 


to  obtain  the  one-way  delay  time  for  the  pulse  between  the  transmitter 
and  the  point  of  contact  with  the  ground. 

For  the  case  of  S(x)  =  1,  these  time  delays  are  plotted  in  Fig. 

10  as  a  function  of  the  take-off  angle  and  parametric  with  the  normalized 
frequency  p,  and  in  Fig.  11  as  function  of  p  and  parametric  with 
The  latter  figure  illustrates  the  caustic  focusing  that  takes  place  at 
minimum  time-delay  and  forms  the  boundary  or  envelope  of  the  sweep- 
frequency  backscatter  echo.  This  envelope,  together  with  that  produced 
by  PQ  =  0,  are  redrawn  in  Fig.  12  for  comparison  with  that  of  Fig.  4, 
which  was  obtained  for  a  parabolic  ion  distribution  by  the  equivalent- 
path  method. 

This  latter  type  of  presentation  of  the  time-delay  results  is 
repeated  for  the  other  conditions  of  the  horizontal  gradient  and  is 
shown  in  Figs.  13  through  l6. 

Finally,  as  a  means  of  predicting  the  approximate  shape  of  the 
intensity  distribution  of  the  backscatter  echo,  a  histographic  plot  is 
made,  from  the  T  vs  curve  (Fig.  10)  of  the  relative  number  of 

rays  that  contribute  towards  each  time-delay  step  or  interval.  Assuming 
an  isotropic  radiation  pattern  for  the  antenna  and,  for  first  order  approxi¬ 
mation,  disregarding  the  effect  of  distance  traveled  by  each  ray,  the  num¬ 
ber  of  rays  (and  hence  the  echo-signal  intensity)  with  any  specified  delay 
interval  such  as  AT  between  T^  and  T^  in  Fig.  IT  will  be  propor¬ 
tional  to  the  angle  enclosed  within  that  interval,  viz.,  (A(3^  +  Apfa) 
covering  both  the  low  and  the  high  rays.  Figure  18  shows  a  sample  histo¬ 
graphic  plot  taken  from  one  of  the  curves  in  Fig.  10  for  p  =  2.25.  The 
time  interval  AT  was  taken  as  0.25  msec.  The  term  "angular  spread" 
refers  to  the  sum  (A£^  +  Aph)  in  degrees. 
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FIG.  10.  VARIATION  OF  PULSE  TIME  DELAY  T  WITH  TAKE-OFF  ANGLE  0o  AT 
VARIOUS  VALUES  OF  NORMALIZED  FREQUENCY  p.  THE  VALUES  OF  T  ARE  OB¬ 
TAINED  FROM  THE  ANALOG  COMPUTER  THROUGH  EQ.  (21).  /3-CHAPMAN  LAYER  OF 
hBax  =  300  ka,  H  =  100  ka,  WITH  NO  HORIZONTAL  GRADIENT  OF  ELECTRON 
DENSITY. 
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FIG.  11.  VARIATION  OF  PULSE  TIME  DELAY  T  WITH  NORMALIZED  OPERATING 
FREQUENCY  p  AT  VARIOUS  VALUES  OF  TAKE-OFF  ANGLE  /3Q  .  Same  layer  condi¬ 
tions  as  those  of  Fig.  10. 
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FIG.  12.  SHAPE  OF  SWEEP- FREQUENCY ,  GROUND- BACKSCATTERED  ECHO  BOUNDED  BY 
ITS  TWO  EXTREME  TIME  DELAYS  AS  FUNCTIONS  OF  p.  /3-Chapman  layer  of  hmax  = 
300  km,  H  =  100  km  with  no  horizontal  gradient.  S(x)  =  1.  Horizontal  lines  within 
the  SFE  are  time  delays  of  rays  meeting  irregularities  at  E-Layer  heights,  with  the 
value  of  ,8  marked  on  each  point.  Data  are  taken  from  ray  tracings  of  Fig.  5. 
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FIG.  13.  SHAPE  OF  SWEEP- FREQUENCY ,  GROUND- BACKSCATTERED  ECHO  BOUNDED  BY 
ITS  TWO  EXTREME  TIME  DELAYS  AS  FUNCTIONS  OF  p.  /3-Chapman  layer  of  h  = 
300  kra,  H  =  100  km  with  horizontal  gradient,  S(x)  =  cob  f(x-b)/R]  with  b/R  =  0  . 
Time  delays  correspond  to  the  ray  tracings  of  Fig.  6. 


FIG.  14.  SHAPE  OF  SWEEP- FREQUENCY ,  GROUND- BACKSCATTERED  ECHO  BOUNDED  BY 
ITS  TWO  EXTREME  TIME  DELAYS  AS  FUNCTIONS  OF  p.  /3-Chapman  layer  of  hm<x  = 
300  km,  H  =  100  km  with  horitontal  gradient,  S(x)  =  cos  [(x-b)/R]  with  b/R  =  +30°. 
Time  delays  correspond  to  the  ray  tracings  of  Fig.  7. 
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FIG.  15.  SHAPE  OF  SWEEP- FREQUENCY ,  GROUND- BACKSCATTERED  ECHO  BOUNDED  BY 
ITS  TWO  EXTREME  TIME  DELAYS  AS  FUNCTIONS  OF  p.  /Z- Chapman  layer  of  hmax  = 

300  km,  H  =  100  km  with  horizontal  gradient,  S(x)  =  cos  [<x-b)/R]  with  b/R  -  +60°. 
Time  delays  correspond  to  the  ray  tracings  of  Fig.  8. 


FIG.  16.  SHAPE  OF  SWEEP- FREQUENCY ,  GBOUND- BACKSCATTEHED  ECHO  BOUNDED  BY 
ITS  TWO  EXTREME  TIME  DELAYS  AS  FUNCTIONS  OF  p.  /3-Chapman  layer  of  h|nax  = 
300  km,  H  =  100  km  with  horiiontal  gradient,  S(x)  =  cos  [(x-b)/R]  with  b/R  =  -60°. 
Time  delays  correspond  to  the  ray  tracings  of  Fig.  9. 
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In  practice,  the  intensity  distribution  will  be  influenced  by 
the  antenna's  vertical  pattern,  resulting  in  the  convolution  of  the  two 
distributions.  For  low-angle  propagation,  such  as  that  employed  in 
the  present  work,  the  effect  will  be  enhancement  of  the  trailing  part 
of  the  echo. 

Further  enhancement  will  be  effected  by  the  finiteness  of  the 
pulse  duration.  This  can  be  seen  by  referring  again  to  Fig.  17,  where 
the  dotted  curve  is  drawn  with  constant  separation  equal  to  the  pulse 
width  t ,  from  the  original  T  -  Pq  curve.  In  this  case  the  angular 
spread  for  a  AT  interval  will  be  that  indicated  by  A'P^  and  A'f^, 
being  the  total  projection  of  the  two  curves  on  the  axis  for  that 

interval . 

E.  ANALYSIS  OF  FACTORS  AFFECTING  SHAPE  OF  BACKSCATTER  ECHO 

In  this  section  the  effects  of  the  various  parameters  of  the  transmis¬ 
sion  circuit  on  the  shape  of  the  ground -backs catter  echo  are  analyzed  to 
gain  insight  into  the  mechanics  of  the  various  shapes  commonly  observed. 
Attention  is  limited  to  the  SFE  bounded  by  the  two  extreme  rays  of  the 
minimum-time  delay  and  of  zero  angle.  The  latter,  which  forms  the 
upper  boundary,  is  not  insignificant  because  of  the  effect  of  the  antenna 
pattern,  particularly  at  higher  frequencies,  as  mentioned  earlier. 

It  is  seen  in  all  the  curves  shown  (Figs.  Il-l6)  that,  because  of 
the  curvature  of  the  earth  and  the  ionosphere,  the  minimum  time  delay, 
particularly  in  the  region  of  the  MUF,  no  longer  bears  the  simple  linear 
relationship  with  frequency  as  under  flat-earth  conditions  (see  Eq.  22 
below) .  This  is  confirmed  in  many  of  the  records  taken  with  the  sounder 
(Sec.  III). 

The  ionospheric  model  assumed  has  a  great  bearing  on  the  shape  of  the 
SFE.  In  the  two  electron-height  distributions  taken- -parabolic  and  the 
p -Chapman --the  difference  is  most  noticeable  in  the  upper  boundary  corre¬ 
sponding  to  low-angle  rays.  Compare  Figs.  4  and  12  and  observe  the 
flatness  of  the  former  as  compared  with  that  of  the  p -Chapman  layer. 

This  behavior  is  to  be  expected,  as  the  two  types  of  distribution  are 
quite  different  at  the  bottom  of  the  layer.  The  constants  for  the  para¬ 
bolic  distribution  used  in  the  determination  of  Fig.  19  were  admittedly 
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FIG.  19.  SHAPE  OF  THE  GROUND- BACKSCATTERED  ECHO, 

BOUNDED  BY  ITS  TWO  EXTREME  TIME  DELAYS  (CORRESPOND¬ 
ING  TO  MINIMUM  TIME  DELAY  AND  TANGENT  RAYS)  AS 
FUNCTIONS  OF  THE  NORMALIZED  FREQUENCY  p.  PARABOLIC 
LAYER  h0  =  140  km,  ym  =  160  km. 

rather  nonrealistic,  but  were  chosen  for  a  best  fit  to  the  Chapman 

distribution  adopted  for  the  ray  tracing.  See  the  three  distributions 

in  Fig.  20  drawn  for  the  same  maximum  ion  density,  i.e.,  10  electrons/ 

3 

cm  at  300 -km  height. 

Comparing  the  SFE  patterns  resulting  from  the  two  parabolic  distri¬ 
butions,  the  effect  of  the  layer  thickness  is,  as  to  be  expected,  to 
increase  the  delay  spread  of  the  SFE.  Hence  a  thin  layer,  such  as  the 

E  or  E  ,  will  return  a  narrow  sweep-frequency  echo.  Similarly,  any 
s 

stratification  or  splitting  in  the  ion  distribution  of  the  layer  would 
be  expected  to  appear  as  separate,  narrow  sweep-frequency  echo.  See  for 
example  the  striations  in  Fig.  21  A-2  at  21. ^40,  Fig.  21  B-2  at  04.20, 
and  Fig.  21  B-U  at  02.20  CMT. 

The  layer  thickness  also  affects,  to  a  certain  extent,  the  slope  of 
the  ainimum-time-delay  part  of  the  echo  as  a  function  of  P  as  can  be 
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FIG.  20.  PLOTS  OF  ELECTRON  DENSITY  VERTICAL 
DISTRIBUTIONS  USED  IN  THE  CALCULATIONS. 


DENSITY  SCALE  IS  RELATIVE  TO  106  ELECTRONS/cc 
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FIG.  21 


FIG.  21 


Part  A-l,  <2,  >3.  Eaatward  propagation  fron  Stanford  daring  the  period  22  -  25  Jane  1962.  Tina,  in  G.M.T. 
progreaaea  fron  right  to  left. 


16  AUGUST  1962 


FIG.  21  B- 2 


19  AUGUST  1962 


G.  21.  SAMPLES  OF  CONTINUOUS- SWEEP- FREQUENCY  RECORDINGS  SHOWING 
VARIOUS  TYPES  OF  ECHOES. 


noted  from  the  following  Eq.  (22)  [Ref.  4]  for  flat-eaAh  geometry  and 
a  parabolic  distribution  in  the  layer, 


# 

where  h  is  the  virtual  reflection  height  corresponding  to  frequency 
.  — y  '  i  ^  ♦ 

f  /  Vl  +  y  A  at  vertical  incidence.  If  this  value  of  h  ,  as  obtained 
c'  »  nr  o  v 

from  the  parabolic-layer  expression  of  Eq.  (A. 8)  in  Appendix  A,  is  substi¬ 
tuted  in  Eq.  (22)  and  all  the  distances  are  normalized  with  respect  to 
the  layer  height  hQ,  we  obtain  the  following  general  expression  for  the 
normalized  oblique  distance  of  minimum-time  delay,  i 


(23A) 


-  S  P  (23B) 

from  which  the  aforementioned  slope  S  can  be  evaluated  in  terms  of  a 

normalized  layer  thickness  y  /h  .  It  will  be  seen  that  the  variation 

m  o 

of  this  slope  S  with  the  parameter  y _/h  is  relatively  slow  (S  varies 

from  2.00  to  4.59  as  (ym/hQ)  increases  from  zero  to  1.0).  Thus,  when 
one  considers  any  particular  SFE  within  a  given  frequency  range,  the 
slope  of  the  leading  edge  of  the  pattern  is  seen  to  be  governed  princi¬ 
pally  by  the  prevailing  critical  frequency  fc  of  the  propagation¬ 
supporting  layer  and  the  height  hQ  of  that  layer. 

Next,  the  effect  of  the  horizontal  variation  of  the  electron  density 
is  to  be  considered.  As  mentioned  earlier,  a  cosine  function  was  chosen 
to  take  account  of  the  effect  of  the  sun's  varying  position  in  the  sky, 
relative  to  the  ray  along  the  path.  Figures  13  through  l6  show  the 
envelope  of  the  SFE  at  different  positions  of  the  subsolar  point  along 
the  path.  Comparison  of  Fig.  13,  for  the  condition  of  noon  at  the 
transmitter,  with  that  of  Fig.  12,  indicating  absence  of  horizontal 
gradients,  shows  that  the  ray  corresponding  to  minimum-time  delay  is, 
again,  virtually  unaffected  by  the  gradient  over  the  path  involved.  The 
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presence  of  a  negative  electron-density  gradient,  however,  causes  some 
of  the  low-angle  rays  to  escape  reflection  at  the  ground,  producing  the 
well-known  layer  tilt  and  making  an  overall  shorter  echo  duration.  This 
effect  is  further  augmented  when  the  subsolar  point  is  located  behind  the 
transmitter,  as  shown  in  Fig.  l6,  resulting  in  shorter  frequency  range, 
since  the  critical  frequency  for  the  remaining  higher-angle  rays  will 
be  lower.  When  the  layer  tilt  is  in  the  opposite  direction,  the  subsolar 
point  being  ahead  of  the  transmitter,  the  slope  of  the  SFE  at  the 
side  is  mostly  affected.  In  both  cases,  however,  the  differences  become 
appreciable  at  positions  of  the  sun  close  to  sunset  and  sunrise. 
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III.  SOUNDING  EXPERIMENTS 


The  ground-backscatter  sounding  was  carried  out  by  means  of  a  novel 
type  of  sweep -frequency  ionosonde  suggested  by  0.  G.  Villard,  Jr.,  and 
built  by  Applied  Technology,  Inc.,  around  a  Hamarlund  SP-600  receiver, 
which  acts  as  a  sensitive  receiver  as  well  as  the  frequency-determining 
element.  The  detailed  design  and  circuit  features  are  given  in  an  un¬ 
published  report.  Figure  22  is  a  block  diagram  of  the  transmitting  and 
receiving  portions  of  the  equipment  when  only  the  exciter  part  of  the 
transmitter  is  used,  giving  a  peak  output  of  about  1  kw.  Higher-power 
outputs  (approximately  10  kw)  are  obtained  by  cascading  the  "final 
amplifier"  of  the  exciter  to  a  "power  amplifier."  Both  of  these  ampli¬ 
fiers  are  L-C  tunable  between  14  and  27  Mc/s,  and  the  tuning  is  made  to 
follow  the  master  frequency  automatically  by  means  of  servo-motors.  A 
special  transmit-receive  (T-R)  arrangement  is  used  for  the  employment 
of  a  single  antenna  for  both  transmission  and  reception. 

The  master  frequency  (f  -  f  )  is  derived  from  the  local  oscillator 
of  the  receiver  (continuous  frequency  fQ)  by  a  weak  capacitive  coupling. 
After  passing  through  an  isolation  amplifier,  the  signal  is  mixed  at 
the  oscillator-mixer  with  a  fixed  frequency  (f^)  of  3-955  Mc/s  (being 
equal  to  the  receiver  first  i-f)  generated  by  a  crystal  oscillator. 

The  mixture  component  (fQ  -  f^)  is  then  fed  back  to  the  receiver  at  its 
rf  input  side,  which  is  also  the  input  for  received  signals  coming  through 
the  T-R  switch  and  the  preamplifier,  and  is  tapped  off  the  rf  output 
through  an  adjacent  isolation  amplifier  and  passed  to  the  driver  section 
of  the  transmitter.  Both  sections  of  the  transmitter,  as  well  as  the 
oxcillator-mixer,  are  held  at  cut-off  in  their  normal  states  and  are 
gated  for  operation  by  the  pulser  during  the  transmission  period.  The 
same  gating  operation  turns  off  the  receiving  preamplifier. 

Frequency  variation  is  effected  by  rotating  the  timing  knob  of  the 
receiver.  Continuous  upward  and  downward  frequency  sweeps  were  made 
by  coupling  to  a  reversible  5-rpm  motor  controlled  by  an  automatic 
switching  mechanism  which  is  actuated  by  two  microswitches  attached  to 
appropriate  points  on  the  rotatable  tuning  scale.  The  frequency  sweep 
was  between  15  and  24  Mc/s  during  a  period  of  approximately  5  min.  in 
each  direction. 
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Because  of  the  continuous  sweep  of  the  frequency,  the  receiver 
would  be  slightly  detuned  with  respect  to  the  returning  echo,  and  this 
effect  would  progressively  increase  with  increasing  echo  delay,  making 
multi-hop  reception  more  difficult.  The  purpose  of  the' slow  sweep  was 
to  reduce  this  effect.  In  addition,  the  receiver  bandwidth  was  increased 
at  the  expense  of  slightly  more  interference. 

As  in  most  conventional  sounders,  the  received  SFE  was  monitored  by 
an  A-scope  and  recorded  directly  on  a  moving  photographic  film  from  a 
second  oscilloscope  whose  beam  was  intensity-modulated  by  the  signal. 

Thus  its  time-base  deflection  was  at  right  angles  to  the  direction  of 
film  travel. 

The  width  of  the  transmitted  pulse  used  was  1  msec.  The  records 
do  not  show  this  directly,  as  it  is  masked  by  the  3-®sec  gate  impressed 
on  the  receiver.  A  delay  of  0.75  msec  takes  place  before  the  generation 
of  the  rf  pulse.  The  PRF  was  about  10  pps. 

At  10-min  intervals,  a  narrow  train  of  1500-km  (l  msec)  range  marks 
is  superimposed  on  the  received  signal.  Both  the  transmitter  pulses 
and  this  train  are  triggered  from  the  same  source. 

The  antenna  system,  which  is  used  for  both  transmission  and  reception, 
is  a  broadside  array  of  eight  rhombics,  which  can  be  made  to  fire  in 
either  the  east  or  the  west  direction  by  appropriate  switching.  Full 
details  of  the  antennas  are  to  be  found  in  Refs.  8  and  9,  but  the  follow¬ 
ing  is  a  summary  of  the  main  characteristics: 


Height  above  ground 
Leg  length 
Apex  angle 

Horizontal  beam  width 
Vertical  beam  width 
Elevation  angle  of  main  lobe 

Some  of  the  records  taken  are 


18.9  meters 
112  meters 
75  deg 
2  deg 
8  deg 

from  8.2  deg  at  25.6  Mc/s  to  l4  deg 
at  15.0  Mc/s 

shown  in  Fig.  21. 
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IV.  DISCUSSION  OF  DATA 


A.  GENERAL 

The  sweep-frequency  ground -backs catter  records,  continuously  obtained 
for  the  frequency  range  of  15  -  25  Mc/s,  show,  on  the  whole,  the  character 
istic  shape  of  the  SFE  derived  by  means  of  ray  tracing  in  Chapter  II. 

This  shape  applies  both  to  the  first  and  to  the  multi-hop  echoes. 

Note  that  the  presentation  of  the  SFE  as  bounded  by  the  two  extreme 
rays  does  not  by  itself  provide  information  as  to  the  relative  intensity 
of  the  various  echo  returns,  which,  as  seen  from  Fig.  l8,  diminishes 
rapidly  at  the  trailing  portion  of  the  echo.  Thus,  details  of  this 
portion  are  sometimes  lost  and  the  upper  boundary  of  the  SFE  will  not 
be  seen.  The  use  of  a  logarithmic  detector  in  the  receiver  helped  a 
great  deal  in  bringing  out  these  details. 

The  presence  of  irregularities  and  departure  from  the  assumed  electron 
density  profile  present  additional  features  to  the  echo  pattern,  appear¬ 
ing  mostly  as  striations  along  the  sweep.  These  are  intensifications  of 
certain  portions  of  the  SFE  at  apparently  random  time  delays.  Most  of 
the  records  appearing  in  Fig.  21  were  chosen  to  show  these  striations. 

“Two  main  types  of  steady,  regular  striations  can  be  identified  in 
relation  to  the  variation  of  their  mean  time  delay  with  frequency: 
l)  frequency-dependent  striations  and  2)  frequency- independent  stria¬ 
tions.  Examples  of  type  l)  are  shown  in  Fig.  21  B-2  around  03.00  GMT 
and  of  type  2)  in  Figs.  21  A-3  at  20. 40,  21  B-l  at  21.50  and  21  B-3  at 
01.50  GMT. 

The  frequency-dependent  striations  are  mostly  of  the  regular  type 
which  follow  a  positive  variation  with  frequency.  There  are  some,  how¬ 
ever,  such  as  those  shown  in  Figs.  21  A-2  at  21. ho,  21  B-l  at  22. *40  and 
21  B-2  at  04.30  CMC  which  exhibit  a  slight  negative  slope.  For  reasons 
that  will  be  explained  in  the  next  section,  this  category  of  echoes 
might  be  more  closely  associated  with  type  2),  which  shows  clear  fre¬ 
quency  independence. 

The  distinction  between  types  l)  and  2)  becomes  less  obvious  when 
the  striation  appears  at  the  trailing  edge  of  the  SFE,  since  this  portion 
of  the  main  echo  is  formed  by  the  low-angle -ray  returns,  which  are 
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themselves  less  frequency-dependent  in  their  time  delay  than  those  at 
the  opposite,  leading  edge,  corresponding  to  minimum-time  delay  (see, 
for  example,  Fig.  4).  As  mentioned  before,  this  is  the  result  of  the 
limiting  effect  caused  by  the  curvature  of  the  earth  and  the  ionosphere. 
Examples  of  these  limit  echoes  are  seen  in  Table  1. 


TABLE  1.  EXAMPLES  OF  LIMIT  ECHOES 


Ea 

St 

West 

Limit  Echoes 

Hours 

Limit  Echoes 

Hours 

22  June 

2300 

17  August 

0030 

23  June 

0230 

17  August 

0430 

23  June 

2210 

18  August 

0200 

24  June 

2030 

19  August 

0150 

The  intensification  of  the  SFE  at  various  slant  ranges  corresponding 
to  their  equivalent -time  delays  can  be  the  result  of  either  l)  time-delay 
focusing;  or  2)  improvement  in  the  reflection  coefficient  at  the  points 
of  incidence  of  the  ray  along  its  path. 

Both  vertical  and  horizontal  gradients  in  ion-density  distribution 
would  result,  at  suitable  angles  and  frequencies  and  with  correct  signs, 
in  some  degree  of  time-delay  focusing.  The  former  gives  rise  to  the 
minimum  time-delay  focusing  appropriate  to  a  layer  with  given  thickness 
and  ionization.  The  same  mechanism  can  be  assumed  to  be  responsible 
for  multiple  intensification  beyond  the  leading  edge  of  the  SFE  if  the 
propagation- supporting  ionosphere  happens  to  be  stratified  or  consisting 
of  small  layers,  separated  perhaps  by  small  ledges  along  the  profile, 
each  contributing  its  own  minimum  time  delay.  The  variation  with  fre¬ 
quency  of  each  of  them  would  be  as  discussed  in  Sec.  II. A  and  governed 
by  Eq.  (22).  The  horizontal  gradients,  on  the  other  hand,  will  result 
in  corrugation  or  ripple  in  the  effective  reflecting  surface  of  the 
ionosphere,  which,  as  postulated  by  Tveten  [Ref.  10],  Hofmann  [Ref.  11] 
and  others,  can  produce  a  certain  degree  of  focusing  in  the  rays. 

Reflection  enhancement  can  take  place  along  the  ray  path  by  the 
presence  of  highly  ionized  irregularities  in  the  F  or  E  layer,  or 
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large  objects  on  the  ground,  such  as  mountain  ranges  [Refs.  12,  13]. 

Such  irregularities  undoubtedly  produce  echoes  that  vary  very  little  with 
frequency,  and  their  relevance  to  the  experimental  results  is  discussed 
in  more  detail  in  the  next  section. 

B.  CONSTANT-RANGE  ECHOES 

1.  General 

This  type  of  echo  has  been  observed  in  this  laboratory  for  some 
time  [Ref.  4],  These  echoes  give  rise  to  frequency-independent  striations 
(type  2),  which,  as  discussed  in  the  previous  section,  represent  echo 
intensification  presumably  resulting  from  irregularities  of  finite  dimen¬ 
sions  located  in  the  ray  path. 

These  striations  are  observed  in  association  with  an  F-layer- 
propagated  echo.  Figure  23  shows  a  plot  of  the  time  and  duration  of 
occurrence  vs  time  delays  of  such  echoes,  extracted  from  the  continuous 
recordings  of  oblique  soundings  in  the  eastward  and  westward  directions 
during  June  and  August,  1962,  respectively.  It  is  seen  that  the  life¬ 
time  of  these  echoes  is  between  a  few  minutes  and  an  hour  or  two. 
Longer-duration  echoes  can  also  be  seen  in  the  plot  but  are,  in  reality, 
trains  of  shorter-duration,  usually  occurring  at  close  intervals  at 
approximately  the  same  slant  ranges.  It  will  thus  not  be  accurate  to 
to  deduce  a  statistical  mean  of  echo  lifetime  based  on  Fig.  23. 

No  detectable  motion  in  these  echoes  is  observed  in  the  records, 
especially  motion  that  would  result  in  a  change  of  the  striation  position 
with  respect  to  the  SFE. 

2.  Ground  Objects  and  F- Layer  Irregularities 

Examination  of  the  slant  ranges  involved  in  the  present  F-layer- 
propagated  backscatter  does  not  reveal  any  definite  or  consistent  ground 
irregularity  that  might  be  taken  as  a  source  of  constant -range  reflection. 
The  possibility  does  exist,  however,  that  the  cluster  of  echoes  observed 
around  3000  km  in  the  eastward  direction  might  be  emanating  from  the 
Florida  coastline,  while  the  Pacific  islands  may  be  responsible  for  the 
more  numerous  western  echoes,  particularly  from  the  Hawaii  region  (4000  - 
5000  km  in  the  direction  of  the  rays).  However,  it  is  difficult  to 
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OCCURRENCE  OF  CONSTANT- RANGE  ECHOES. 


explain  the  relatively  short  and  sporadic  appearance  of  the  echoes  on 
that  basis. 

Though  of  little  significance  to  the  case  of  oblique  transmissions, 
analysis  of  the  regular  vertical -incidence  ionograms  taken  by  the  C-2 
sounder  at  Stanford  shows,  for  the  same  periods,  absence  of  any  echoes 
that  might  originate  from  some  ion-concentration  irregularities  in  the 
F  layer  overhead. 

The  effect  of  a  reflecting  irregularity  in  the  F  layer  can  be 
followed  in  any  of  the  ray  tracings  presented  above,  although  the  actual 
time  delay  cannot  be  obtained  accurately  since  the  computed  time  delays 
have  only  been  integrated  over  the  full  path  of  the  ray.  (An  approxi¬ 
mation  was  made  for  the  E-layer  backs cattering  by  subtracting  the  delay 
in  the  100-km-high,  non-deviative  portion  from  the  total  delay.  See 
Sec.  IV .A. 3.)  However,  two  features  of  the  echoes  that  would  be  obtained 
can  be  deduced  from  a  study  of  the  ray  tracings: 

1.  The  time  delay  of  such  echoes  would  be  fairly  independent  of 
frequency,  though  they  would  exhibit  a  slight  negative  slope, 
which  would  generally  be  less  than  that  for  the  E-layer  back- 
scattering. 

2.  The  constant-range  echoes  thus  produced  would  invariably  extend 
well  beyond  the  edge  formed  by  the  minimum-time -delay  portion  of 
the  SFE,  since  some  reflection  from  the  irregularity  can  conceiv¬ 
ably  occur  within  the  skip  zone.  The  amount  of  extension  would 
depend  on  the  relative  location  of  the  irregularity  in  the  F  layer 
and  on  its  ionization  intensity.  Such  appreciable  protrusions 
are  not  observed  in  the  recordings. 

The  foregoing  discussion  regarding  the  F  layer  is  based  on  the 
assumption  that  the  reflecting  irregularities  are  of  such  high  local 
electron  densities  that  they  will  possess  sharp  gradients  and  thus  produce 
forward  scattering.  F-layer  irregularities  of  milder  gradients  would 
allow  at  least  partial  penetration  of  the  rays,  thereby  suffering  re  - 
tardation,  and  would  therefore  exhibit  frequency  dependence  in  their 
time  delay. 

The  vertical-incidence  ionograms  which  were  examined  also  show 
approximately  equal  sporadic-E  activity  overhead  during  the  two  periods, 
yet  the  backscatter  sounding  shows  definite  preponderance  of  E-layer 
echoes  at  about  1000  km  when  the  observation  was  directed  inland 
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(Fig.  23.A).  This  increase  may  well  be  a  reflection  from  the  Rocky 
Mountain  range  [Ref.  15].  On  the  other  hand,  the  following  may  be  an 
alternative  cause  of  these  close-range  echoes. 

3.  Irregularities  in  the  E  Layer 

Peterson  et  al  [Ref.  16]  and  Gallagher  [Ref.  17]  observed  field- 
aligned  irregularities  in  the  E  region  (centered  around  the  geomagnetic 
north  from  Stanford),  which  produced  strong,  but  rapidly  fading,  direct- 
backscatter  echoes  resulting  from  the  incidence  of  the  radar  wave  on 
the  elongated  ionization  at  right  angles  to  its  long  axis  of  symmetry. 
These  middle-latitude  irregularities  may  be  associated  with  extensions 
of  the  polar  or  equatorial  type  of  sporadic-E  to  middle  latitudes. 

The  slant  range  at  which  this  pnenomenon  has  been  observed  around  Stan¬ 
ford,  in  the  absence  of  ionospheric  bending,  depends  on  the  direction 
of  propagation,  and  ranges  from  about  200  km  due  magnetic  north  (declin¬ 
ation  17.5°E)  to  more  than  1000  km  on  either  side  of  this  direction. 

This  phenomenon  certainly  cannot  account  for  all  of  the  constant-range 
echoes  observed  in  the  present  experiments,  some  of  which  extend  beyond 
6000  km,  but  one  can  picture  a  ray  meeting  the  axis  of  the  irregularity 
orthogonally  at  the  E-layer  level  after  being  reflected  by  the  upper 
F  layer.  With  this  geometry,  and  Stanford  as  center,  it  is  possible  to 
meet  the  geomagnetic-field  lines  at  right  angles  for  distances  extending 
beyond  8000  km.  It  is  shown  in  Appendix  D  that  normal  intersection  with 
a  field  line  F  of  positive  dip  angle  I  takes  place  when  the  vertical  or 
elevation  angle  (with  respect  to  the  horizontal)  of  the  downcoming 

ray  satisfies  the  following  simple  relation: 

tan  S  =  cos  (l80  -  a)  cot  I,  (23) 

3? 

with  a  being  the  angle  between  the  two  meridian  planes  containing  the 
ray  and  the  geomagnetic  field  line  at  the  point  of  intersection  P . 

If  the  directions  of  these  planes  relative  to  the  magnetic  north  are 
A  and  D  (declination)  respectively,  the  angle  a  will  be: 

for  westward  propagation,  a  -  180  -  (A  -  D).  (24) 

for  eastward  propagation,  a  =  180  -  (A  +  D).  (25) 
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With  values  of  I  and  a  corresponding  to  points  P  along  the 
two  directions  of  propagation,  a  plot  is  presented  in  Fig.  24  giving  the 
necessary  value  of  p^  at  these  points  to  satisfy  the  perpendicularity 
condition.  The  negative  values  of  p^  to  the  east  indicate  that  the 
ray  meets  the  line  of  force  in  its  upward  direction. 


It  will  be  noticed  from  Fig.  24  that,  under  this  assumption  of 

meeting  the  irregularities  in  a  downcoming  direction  of  the  ray,  the 

westward  propagation  requires  appreciably  larger  values  of  the  angle 

Pp  to  achieve  perpendicular  intersection.  Inspection  of  the  rays  traced 

for  a  layer  in  the  absence  of  a  horizontal  gradient  (Fig.  5)  shows  that 

the  angle  of  the  downcoming  ray  meeting  an  irregularity  at  a  height  of 

100  km  at  a  given  range  does  not  vary  materially  from  10  deg.  The 

horizontal  lines  shown  embedded  in  the  SFE  pattern  of  Fig.  12  are  one- 

way  time  delays  calculated  from  these  tracings  between  the  transmitter 
* 

Calculated  from  the  data  provided  by  the  geomagnetic  maps  in  [Ref,  18], 
Values  at  E-layer  height  are  considered  to  be  approximately  the  same 
as  those  given  for  ground  level. 
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and  points  P  distant  1750,  2000,  and  2250  km,  respectively,  from  it 
and  at  100-km  height.  The  values  of  shown  are  average  gp  of 

all  the  rays  that  strike  each  particular  point  P  at  that  height. 

Thus,  it  is  seen  that,  although  the  time  delays  of  the  rays 
striking  the  irregularity  at  the  point  P  remain  fairly  constant, 
particularly  at  the  more  distant  points,  the  average  value  of  the  required 
angle  g  appears  to  he  too  high  for  the  easterly  direction  of  propaga¬ 
tion  and  too  low  for  the  opposite  direction. 

Possibly  one  way  of  reconciling  these  discrepancies  is  to  consider 
the  time  of  appearance  of  these  constant-range  echoes  on  a  local  basis, 
as  shown  in  the  histograms  of  Fig.  25.  It  is  found  that  the  mean  of  this 
time  is  about  the  same  in  either  direction,  with  a  tendency  towards  the 
early  afternoon  for  the  westward  propagation  and  towards  the  late  after¬ 
noon  for  the  eastward  propagation.  This  approximately  constant  mean  time 
of  appearance  in  the  afternoon  indicates  that  the  horizontal  gradients 
of  the  ion  densities  in  the  F  layer  are  of  opposite  signs  in  the  two 
directions,  resulting  in  layer  tilts,  the  senses  of  which  do  tally  with 
the  inequality  of  g^  for  given  ranges  in  these  directions  (Fig.  24). 
However,  the  layer  tilts  produced  by  the  cosine  gradient  assumed  for 
the  SFE  of  Fig.  15,  which  corresponds  to  the  westerly  direction,  does 
not  seem  to  bring  about  a  substantial  increase  in  the  angle  g^  (the 
mean  values  are  also  shown  for  each  ground  range). 

There  are  two  features  of  the  constant -range  echoes  that  the  ray 
tracings  have  brought  out,  under  the  above  assumption  of  the  E-layer 
origin  of  the  backscattering  irregularities.  These  are  l)  their 
slight  negative  slope,  diminishing  rapidly  at  longer  ranges;  and  2) 
their  extension  slightly  beyond  the  minimum-time-delay  boundary  (by 
an  amount  equivalent,  for  any  given  frequency,  to  the  oblique  distance 
of  the  irregularity  from  the  ground-reflection  point  which  the  inter¬ 
cepted  ray  would  have  reached).  Both  of  these  characteristics  are 
observed  in  varying  degrees  in  many  of  the  echoes  found  in  the  recordings, 
as  can  be  seen,  for  example,  in  Figs.  21  B-2  at  04.50  and  21  B-4  at  01.50 
GMT,  though  many  others  do  not  exhibit  these  features  but  still  appear 
frequency  independent. 
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a.  Eastward  Propagation. 


itesi  Tint 


b.  Westward  Propagation. 


FIG.  25.  STATISTICAL  PLOT  OF  THE  FREQUENCY  OF  OCCURRENCE  OF  CONSTANT 
RANGE  ECHOES  AGAINST  TIME  OF  OCCURRENCE  PREVAILING  LOCALLY  AT  THE 
POSITION  OF  THE  IRREGULARITY. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  ray  analysis  has  shown  that  the  shape  of  the  SFE  pattern  (delay 
time  vs  frequency)  derived  from  either  of  the  two  methods  discussed  in 
Sections  II. C  and  II. D,  respectively,  agrees  well  with  that  observed  in 
normal  sweep-frequency  sounding.  From  a  study  of  the  shape,  certain 
properties  of  the  propagation-supporting  ionosphere  can  be  surmised, 
such  as  layer  thickness  and  the  presence  of  extensive  horizontal  grad¬ 
ients  . 

It  has  further  been  shown,  and  graphically  demonstrated  by  an 
extension  of  Peterson's  construction  to  curved-earth  geometry,  that  the 
tangent  ray  makes  a  definite  contribution  at  the  trailing  edge  of  the 
above  pattern,  and  varies  least  with  frequency,  approaching  frequency 
independence  between  wide  limits.  This  is  noteworthy  in  respect  to 
the  constant-range  striations  with  which  this  study  has  been  chiefly 
concerned.  The  study  of  the  two  types  of  striations  appearing  within 
the  SFE,  classified  with  respect  to  the  variation  of  their  time  delay 
with  frequency,  has  revealed  that  those  which  vary  positively  are  the 
result  of  echo-delay  focusing  produced  by  suitable  vertical  or  horizontal 
gradients  in  the  ionosphere.  On  the  other  hand,  highly  reflective 
irregularities  (such  as  sharp  gradients)  whether  midway  along  the  path 
of  the  ray  (in  the  F  region)  or  at  its  far  end  (in  the  E  region  or  the 
ground)  can  produce  echoes  whose  delay  times  are  fairly  independent  of 
the  frequency. 

As  far  as  the  latter  type  of  echo  is  concerned,  analysis  of  the 
backscatter  data  collected  for  the  two  directions  of  propagation  has 
shown  that,  of  the  three  possibilities,  the  E-layer  origin  of  the  ex¬ 
pected  irregularities  is  the  most  likely.  It  is  also  quite  feasible 
that  these  irregularities  are  in  line  with  the  geomagnetic  field,  al¬ 
though,  for  the  layer  model  adopted,  some  kind  of  layer  tilt  must  be 
invoked  if  the  maximum- reflection  condition  is  to  be  fulfilled.  Their 
predominantly  daytime  appearance  seems,  a  priori,  to  dissociate  such 
echoes  from  the  nocturnal  type  observed  by  Gallagher  [Ref.  17]  in  the 
same  layer,  but  this  finding  is  far  from  being  conclusive  because  of 
the  short  period  of  observation  of  the  constant-range  echoes  in  the 
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present  study,  during  a  great  part  of  which  night  propagation  via  the 
F  layer  was  absent. 

Frequency-sweep  sounding  has  brought  to  light  certain  types  of 
irregularities  that  would  not  be  distinguishable  in  normal  fixed- 
frequency  backscatter  sounding,  nor  ordinarily  in  multifrequency  opera¬ 
tion.  Further  investigation  needs  to  be  carried  out  to  establish  their 
morphology  and  azimuthal  extent,  and  to  find  out  more  about  their  loca¬ 
tion  and  apparent  motion. 
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APPENDIX  A.  DERIVATION  OP  THE  EQUIVALENT  GROUP  PATH  AS  A  FUNCTION  OP 
THE  TAKE-OFF  ANGLE  FOR  CURVED  GEOMETRY 


1.  In  Terms  of  Virtual  Height 

Fig.  1,  the  sine  law  for  the  triangle  QAB  gives : 


P'  =  2(R  +  h • ) 
=  2(R  +  h') 


sin  6 
cos  6 

o 

cos  (8  +  0  ) 
O  0 

cos  8 
Ko 


(A.l) 


and  since 


cos 


(8  +  0  )  =  cos  8  cos  0 
o  o  o  o 


sin  8  sin  0  , 
o  o 


then 

P1  =  2(R  +  h')  (cos  0q  -  tan  8Q  sin  0^).  (A. 2) 

Again,  from  the  sine  law  for  the  triangle  QAB,  and  letting  r  = 

R/(R  +  h'),  we  get 


and 


sin  0  =  r  cos  8 
o  'o 


tan  8  0  =  r  sin  8 

'o  o  o 


Also 


or 


cos  0  =  (l  -  sin2  0 
*o  ro 

2  2  a  vl/2 

=  (1  -  r  cos  8  )  • 

o' 

From  (A. 2),  (A. 4)  and  (A. 5),  we  get 

P'  =  2(R  +  h')  [(1  -  r2  cos2  0O)1//2  -  r  sin  p  ] 

P'  =  2  { [(R  +  h')2  -  R2  cos2  80]^2  -  R  sin  0Q} . 


(A. 3) 


(A.  4) 


(A. 5) 


(A.6) 

(A. 7) 
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2.  In  Terms  of  Layer  Parameters 

If  a  parobolic  distribution  of  electron  density  is  assumed  for 
the  layer,  then  for  vertical  incidence  of  a  wave  at  frequency  f 


h'  =  h  + 
o 


,  fl  +  (f  /f  )  “) 

2  ym  in[l  -  (//mJ  (A,8) 


where:  hQ  =  height  of  the  bottom  of  the  layer 

y  =  half  thickness  of  the  layer 
m 

f  =  plasma  frequency  or  critical  frequency  at  vertical  incidence 
At  oblique  incidence,  and  when  the  effect  of  the  curvature  of 
the  ionosphere  at  the  point  of  "reflection"  of  the  wave  of  frequency 
f  ^  is  neglected, 


f  =  f  ,  cos  0  . 
v  ob  o 


(A. 9) 


o  ■  <  W 


__  =  p  COS  0 

f  o 

c 


(A. 10) 


and,  from  Eq.  (A. 5), 


Y-  =  P  [1  -  r2  cos2  PQ]1/2. 
c 


(A. 11) 


When  the  above  is  put  in  Eq.  (A. 8)  and  is  expressed  in  terms  of 
(R  +  h'): 

1  2  2  ,1  1  +  P(1  -  r2  cos2  p  )2 

(R  +  h')  =  (R  +  hQ)  +  |  ym  p(l  -  r  cos^  ^  in - - - r-^r 

1  -  p(l  -  r  cos ,  PQ)2 


(A. 12) 

Remembering  that  r  is  a  function  of  (R  +  h ' )>  we  evaluate  the 
latter  by  successive  approximations -using  a  digital  computer. 
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3.  Rate  of  Variation 


Putting  Ri  =  (R  +  h1)  in  Eq.  (A. 7)  and.  differentiating  with 
respect  to  we  obtain 

o  p  p  -i  2 

[£(R7  -  R  cos  p  )"2  (2R,  +  2R  sin  p„  cos  6  )  -  R  cos  B  ] 

I  0  1  ap  0  0  o 


dP’ 


2  2  2  2“ 

-TTri  +  R  sin  B  cos  p  )  (R  -  R  cos  p  )  2  -  R  cos  p 
dp  Ko  o'  i  Ko'  Ko- 

h  o 


(A. 13) 


For  tangential  propagation  PQ  *»  0,  hence 


dP> 

W 

o 


-  2 


[Ri  w0  (Ri  -  r2)4  -  R] 

“H  ffi 

J  d^0  ’ 


f(R  -  R)(R  +  R)*]4  dR 

'  2L - 2 - J 


(A. 14) 


Eq.  (A.l4)  can  further  be  approximated  by  putting  (Ri  +  R)  *  2R;  we 
then  obtain 


dP’  ^  /  2R  \  £  dh^_ 

dp  *  \  h '  )  dp 
Ko  o 


2R 


(A. 15) 


At  low  values  of  p  the  virtual  height  h'  tends  to  be  con- 
o 

stant,  and  the  first  term  of  (A. 15)  becomes  negligible.  Therefore 
dP'/dp  -»  (-  2R)  as  p  ->  0. 

'  O  O 
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APPENDIX  B.  GEOMETRY  OF  TRANSMISSION  PATH  WITH  RESPECT  TO  SUN'S 
ZENITH  ANGLE 


Q 


In  the  above  sketch,  T  represents  the  transmitter  and  Q  the  pro¬ 
jection  on  the  earth's  surface  of  a  point  q  in  the  ray  path.  S  is 
the  subsolar  point,  and  b,  x,  and  y  are  distances  on  the  surface  of 
the  earth  (of  radius  R) .  The  sun's  zenith  angle  is  therefore  given  by: 

X  =  (y/R)  (B.la) 

and 

cos  X  =  co 6  (y/R)  (B.lb) 

From  spherical  trigonometry 

cos  ^  =  cos  g  cos  ^  +  sin  ^  sin  ^  cos  y  (B.2) 

Two  limiting  cases  are  considered: 

1.  The  subsolar  point  along  the  x  direction,  i.e.,  in  the  plane  of 
the  ray  path,  7=0, 

cos  (y/R)  =  cos  (x/r)  cos  (b/R)  +  sin  (x/R)  sin  (b/R) 

cos  X  =  cos  [(x/R)  -  (b/R)]  (B.3) 

2.  The  subsolar  point  in  a  meridian  plane  orthogonal  to  that  of  the 
ray  path  y  =  n/2 

cos  X  =  cos  (x/R)  •  cos  (b/R)  (B.4) 
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APPENDIX  C.  DERIVATION  OF  ECHO  GROUP  TIME  DELAY 


In  the  ray  tracings  of  Fig.  5  et  seq.,  which  are  plotted  in  terms 
of  the  ground  distance  x  (from  the  transmitter)  and  the  true  height 
h,  the  length  of  a  elemental  portion  of  the  traced  curve  is 

ds  =  [(dx)2  +  (dh)2]7  (C.l) 


Hence 


(C.2) 


However,  although  the  tracings  give  the  coordinates  of  any  point  in 
the  ray  path  in  terms  of  x  and  h  correctly,  the  Cartesian  system 
does  not  yield  the  true  shape  or  the  true  length  of  the  ray  path,  which 
follows  the  curvature  of  the  earth.  With  this  approximation  in  mind, 
the  group  delay  is  derived  from  the  ray  equation  as  follows: 


Let  t  =  time  elapse  in  pulse  propagation 

T  =  total  delay  of  pulse  for  one  hop  and  one  way 
v  =  group  velocity  of  wave 

o 

p  =  phase  refractive  index 
p  *  =  group  refractive  index 
c  =  wave  velocity  in  free  space 

Then 

dt  =  ds/v  =  (p'/c)  ds  (C . 3) 

o 


Under  the  ionospheric  conditions  assumed  in  this  work,  the  effect 
of  the  earth's  magnetic  field  and  of  collisional  damping  can  be  neglected, 
so  that: 

VJL  =  1/PL '  (C.4) 

and 

dt  =  (l/cp)  ds 

=  (l/cp)(ds/dx)  dx.  (C.5) 
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-ilm 

path 

dx 


(C.6) 


Prom  Eq.  (l4)  in  the  text, 

h  =  0 


x  =  0L-L  " 


(dh/dx)‘ 


P(h)  •  S(x)/p' 


r 


dx 


(C.7) 


The  terms  in  Eq.  (C . 7)  are  tapped  off  the  ray-tracing  program. 
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APPENDIX  D.  DETERMINATION  OF  THE  ANGLE  OF  THE  DOWNCOMING  RAY  MEETING  A 
GEOMAGNETIC  FIELD  LINE  AT  RIGHT  ANGIES 


In  Fig.  26,  which  is  drawn  for  westward  propagation,  in  the  plane 
of  the  paper,  let  P  be  the  point  of  intersection  of  the  ray  CP  in 
this  meridian  plane  with  the  geomagnetic  field  line  F  contained  in 
another  meridian  plane  making  an  angle  a  with  the  former.  In  the 
present  work,  the  location  of  P  is  at  a  height  of  100  km.  Also,  let 
P  and  I  be  the  angles  which  the  two  vectors  make  with  the  horizontal 

p  NMIBIM  FU«  OF  F  (KMMETIC 


FIG.  26.  GEOMETRY  OF  THE  NORMAL  INTERSECTION  OF  A  DOWNCOMING  RAY 
WITH  A  GEOMAGNETIC  FIELD  LINE  LOOKING  WEST  OF  STANFORD. 

In  a  plane  cutting  the  two  aforementioned  meridians  orthogonally, 
a  triangle  ABC  is  formed  from  its  intersection  with  the  ray  and  the 
magnetic  field,  with  angle  ABC  equal  to  a,  so  that 


AC2  =  AB2  +  BC2  -  2  AB  BC  cos  a. 


(D.l) 


In  the  two  right-angled  triangles  with  P  as  the  common  apex 
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and 


IS  ■  cos  I 


(D.2) 


BC  ■  PC  cos  p. 

Combining  the  above  three  equations 

AC2  =  PA?  cos2  I  +  PC1*  cos2  p  -  2  PA  •  Pi?  cos  a  cos 
*• 

Also  in  triangle  APC 

AC2  =  PA2  +  PC2. 

Equating  (D.4)  and  (D.5)  and  rearranging,  we  get 

PA2  sin2  I  +  PC2  sin2  p  =  -  2  PA  •  PC  cos  a  cos  p 

Again,  in  the  two  right-angled  triangles  having  PB  as 
side, 

PB  =  KT  sin  p  =  SI  sin  I. 

From  (D.6)  and  (D.7)  we  obtain  the  relationship 

ten  p  =  -  cos  a  cot  I. 


(D.3) 


p  cos  I.  (D.4) 


(D.5) 


cos  I.  (D.6) 
the  common 

(D.7) 


(D.8) 
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